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ABSTRACT 
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^ ; We have investigated time evolutions of power spectra of density fluctuations 

for long time after the first appearance of caustics in the expanding one- 
O ' dimensional universe. It is found that when an initial power spectrum is 
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sale-free with a power index n, a self-similarity of the time evolution of the 
power spectrum is achieved. We find that the power spectrum can be separated 
roughly into three regimes according to the shape of the power spectrum: the 
linear regime (k < k n i : the regime 1), the single- caustic regime (k n i < k < k sn i 
: the regime 2), and the multi-caustics regime (k > k sn i : the regime 3). The 
power index of the power spectrum in each regime has the values of n, — 1, and 
\i which depends on n, respectively. Even in the case of an initial power-law 
^ ■ spectrum with a cutoff scale, there might be the possibility of the self-similar 

evolution of the power spectrum after the appearance of the caustics. It is 
found, however, the self-similarity is not achieved in this case. The shape of the 
power spectrum on scales smaller than the cutoff scale can be separated roughly 
in two regimes: the virialized regime (k cut < k < k cs : the regime 4), and the 
smallest-single-caustic regime (k > k cs : the regime 5). The power index of the 
power spectrum is v which may be determined by the distribution of singular 
points in the regime 4. In the regime 5, the value of the power index is —1. 
Moreover we show the general property about the shape of a power spectrum 
with a general initial condition. 

Subject headings: cosmology:theory-large scale structures-selfsimilar 
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1. INTRODUCTION 

Formations of the large scale structures in the expanding universe is one of the most 
important and interesting problems of the cosmology. It is generally believed that these 
structures have been formed owing to gravitational instability. Hence it is very important to 
clarify evolutions of density fluctuations by the gravitational instability. Here we consider 
density fluctuations of collisionless particles as dark matters, our interest being mainly 
concentrated on effects of the self-gravity. 

The density fluctuations are often characterized by their power spectrum (or two-point 
correlation function). Thus, it is very interesting and important to investigate time 
evolutions of the power spectrum. Furthermore, whether the evolutions of the power 
spectrum satisfy the self-similarity or not is especially important. If the self-similarity is 
satisfied, we can describe completely the time evolution of the power spectrum without 
numerical simulations because we can estimate the growth rate of the power spectrum on 
all scales. Many authors have investigated the self-similarity of the two-point correlation 
function (or the power spectrum). Davis & Peebles (1977) investigated the self-similarity of 
the two-point correlation function by using the BBGKY equations. It is well known that 
we can make a similarity transformed BBGKY equations because the gravity is scale-free. 
They showed the existence of the self-similarity by integrating the BBGKY equations 
numerically under the assumptions as shown below. The BBGKY equations have the 
hierarchical structure, that is, the time evolution of the iV-th order correlation function 
includes the (N + l)-th order correlation function. So, the cutoff of the hierarchy is needed 
in order to close these equations. Davis & Peebles assumed that the three-point correlation 
function can be expressed by the product of the two-point correlation functions and that 
the skewness of the velocity field is equal to zero. Thus the existence of the self-similarity 
is still uncertain because we do not know whether these assumptions are correct or not. 
On the other hand, by using iV-body simulations, time evolutions of the power spectrum 
can be calculated directly. Recent works, for example, are Colombi,Bouchet & Hernquist 
(1996); Couchman & Peebles (1997); Jain, Mo, & White (1995) and so on. They showed 
that the self-similar evolution of the power spectrum can be satisfied when the initial power 
spectrum is scale-free. However the relation between the index of the power spectrum and 
the mean relative velocity is still uncertain quantitatively. Furthermore, only the scale-free 
cases are investigated. In this paper we would like to investigate not only the scale-free 
cases but also the cases that the initial power spectrum obeys the power-law with a cutoff. 
Hereafter we call this case cutoff-case. In the cutoff-case, the power spectrum obeys the 
power law even on scales smaller than the cutoff scale after the first appearance of caustics. 
For example, in the one-dimensional system, The value of the power index is —1 which can 
be derived according to the catastrophe theory (Gouda & Nakamura (1988,1989)). Kotok 
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& Shandarin (1988) also studied the nonlinear spectra in the cutoff case and showed that 
the value of the power index is -1. Of course, the power spectrum does not evolve in the 
self-similar form before the first appearance of caustics. However it has not been certain 
whether the time evolution is self-similar or not after the first appearance of caustics. 
Yano & Gouda (1997) showed the possibility of the self-similarity in this case. Therefore 
we investigate the self-similarity not only in the case of the scale-free case but also in the 
cutoff-case. Then we have to investigate in detail the evolutions of the power spectrum 
on scales smaller than the initial cutoff scale and also those on larger scales. Thus we 
need a wide dynamic range in the wave number space. However, we cannot get such a 
wide dynamical range in the numerical simulation of three-dimensional systems due to the 
limit of the resolution. Then, as a first step, we consider one-dimensional sheet systems. 
In the one-dimensional system we can get a wide range of resolution for calculating time 
evolutions of the power spectrum. Furthermore we can have the numerical method for the 
evolution of the power spectrum with a good accuracy in the one- dimensional sheet systems 
as shown in §2. The purpose of this paper is to investigate the time evolution of the power 
spectrum in the one-dimensional systems as a first step. We are especially interested in the 
self-similar evolution of the power spectrum. Therefore, we consider only the background 
universe which does not have a characteristic scale, as a spatial curvature scale. That is, we 
consider the Einstein-de Sitter universe. 

In §2.1, we will briefly show the numerical method. We show in §2.2 and §2.3 the 
numerical results of scale-free case, cutoff-case, respectively. Finally, we will devote to the 
conclusion and discussions in §3. 



We investigate time evolutions of the power spectrum and its self-similarity by using a 
numerical method. In this paper we consider the one-dimensional system that many plane 
parallel sheets move only a perpendicular direction to the surface of these sheets. When two 
sheets cross, they are allowed to pass through freely each other. In this sheet system, there 
is an exact solution until two sheets cross over as follows (Sunyaev & Zel'dovich(1972), 
Doroshkevich et.al.(1973)): 



2. 



NUMERICAL RESULTS 



2.1. 



Numerical method 



x 



v 



q + B 1 (t)S 1 {q) + B 2 (t)S 2 (q), 
B 1 (t)S 1 (q) + B 2 (t)S 2 (q), 



(1) 
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where q and x are the Lagrangian and the Eulerian coordinate, respectively. Here, Si(q), 
and S 2 (q) are arbitrary functions of q. Bi(t), and B 2 {t) are the growing mode and 
the decaying mode of linear perturbation solutions, respectively. Since we consider the 

3 

Einstein-de Sitter universe, -Bi(t) = a and B 2 (t) = a~ 5 , where a is a scale factor of the 
universe. We can compute the crossing time of all neighboring pairs of sheets. We use a 
shortest of these crossing time as a time step. Then, we can compute the new positions 
and velocities for all sheets at this crossing time. After two sheets cross, we exchange the 
velocities of just crossed two sheets. Then we obtain again S\(q), S 2 (q), and therefore exact 
solutions as follows: 

3 2 
Si(q) = -a^(x - q) + ^a~V 

S 2 (q) = -a* (a; - q) - -a^ah. (2) 

These new exact solutions can be used until two sheets cross over. In this way we obtain 
the exact loci of the sheets by coupling these solutions. This is the numerical method with 
good accuracy because we connect the exact solutions. Through this paper, we use 2 13 
sheets for numerical calculation. A period boundary is fixed by a length of 2tc. Therefore, 
the smallest wave number is 1. 



2.2. Scale-free spectrum case 

In Figs. 1(a) and (b), we show the time evolutions of the power spectrum with the 
scale-free initial power spectrum given by 

P(k,t ini )<xk n , (3) 

where n is the power index of the initial spectrum. It is found that the following results are 
satisfied for — 1 < n < 4 . The condition of n > — 1 means the hierarchical clustering in 
the one-dimensional system. We are interested in this case ( n > — 1 ) for "scale- free" case. 
One of the reasons why we consider this case is that the hierarchical clustering picture is 
expected in the real world. Another reason is that the results for n < — 1 is similar to those 
in the case of the single wave case. See also eq.(14). Then the results for n < — 1 can be 
reffered from the results shown in the regime4 and 5 in the "cutoff case" (see §2.3). The 
condition of n < 4 is required because we consider the situation in which the non-linear 
mode coupling from higher k to smaller k can be neglected and so the linear perturbation 
theory can be hold on small k( Peebles (1980), Shandarin & Melott (1990), Gouda (1995)). 
In the following, we show the case of n — 1 and 2 for example. Fig. 1(a) and(b) are the 
n = 1 case and the n = 2 case, respectively. We are interested in the hierarchical clustering 



- 5 - 



picture. Therefore we chose some examples that have the power index larger than — 1. We 
can obtain the same results qualitatively in the cases of the other indexes. Here tj„j is the 
initial time, and the initial scale factor of these two cases are 0.1 and 0.01, respectively. We 
normalize the scale factor as follows; a = 1 when the first caustic has appeared. The phases 
of the initial fourier spectrum are given in random. We have averaged 50 samples. Here, we 
scale the wave number and power spectrum as follows: 

where k n i is defined by 

1 fk n l(t) 

— P(k,t)dk = l. (5) 

Zir Jo 

In the regime of k < k n i, the power spectrum grows according to the solution of linear 
perturbation, that is, the power spectrum satisfies the relation P(k,t) oc a 2 . Therefore, k n \ 
is proportional to a~ 2 ^ n+1 \ P sca ie(t) is defined by 

P scale (t) = P(k nl (t),t). (6) 

Then P SC aie{t) is proportional to a 2 /^ oc a 2 ^ n+1 \ This scaled power spectrum P*(/c*) 
is shown in Fig. 2 (a) for the case of n — 1 . The same one but for the case of n = 2 
is shown in Fig.2(b). We can see the coincidence of each power spectrum at each time 
with good accuracy. This means that the self-similarity is achieved in the scale-free case. 
Furthermore we can see the three different power indexes in three regimes: the linear regime 
(k < k n i : the regime 1), the single-caustic regime (k n i < k < k sn i : the regime 2), and the 
multi-caustics regime (k > k sn i : the regime 3). Fig. 3 shows the schematic general power 
spectrum at a time in the scale-free case. The value of the power index in the linear regime, 
of course, remains unchanged. We explain the power index of the power spectrum in the 
other two regimes: the single- caustic regime, and the multi-caustics regime in the following 
subsections. 



2.2.1. Single- caustic regime (the regime 2) 

In the single-caustic regime (k n i < k < k sn i : the regime 2) we can see that the power 
index value becomes —1. The reason is as follows; The density perturbation with the scale 
which just entered the non-linear regime at a given time, that is, k n i, makes "caustics". 
Strictly speaking, the caustics would be appeared if the initial spectrum was smoothed 
bellow the scales just entered the non-linear regime. Therefore, we must notice that the real 
caustics cannot be observed in the scale-free case because of smearing by the small scale 
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fluctuations. However we call these scale on k n i < k < k snl} "single-caustics regime" in this 
paper. After the first appearance of "caustics", the value of the power index of the power 
spectrum on the scales smaller than that scale is —1 (Gouda & Nakamura (1989), Kotok & 
Shandarin (1988)): Here we briefly show why the value of the power index is —1 after the 
first appearance of caustics. 

The fourier spectrum 5k of the density fluctuation is given by 

S k = J 5{x)e lkx dx. (7) 

Then the density is given by 

= (8) 

I dq I 

where po is the mean density of the universe. At the Lagrangian singular point go, the 
following relation is satisfied. 

(f ),„ = 0. (9) 

At the singular point q , the density p(x) diverges. Therefore, the Eulerian coordinate x 
can be written 

x = x + (3(q - g ) 2 + 0((g - g ) 3 ). (10) 

For simplicity, we can put xo = go = without losing generality. Therefore we can express 
x = (3q 2 around the caustic. The density around the caustic is expressed by 

pO) = Pol^T 1 oc (px)~* (11) 

This density profile determine the proper index of the power spectrum. This type of 
singularity is called A2 type according to the catastrophe theory (Gouda & Nakamura 
(1988,1989)). Eqs.flTUD, and (O) can be generally satisfied around the singular points in 
the multi-stream flow regimes where the Zel'dovich solution cannot hold. The A2 type of 
singularity is the only stable type in the one- dimensional system and their eqs.(|l0|), and ( 11) 



remain in the multi-stream flow regions after the first appearance of caustics. Here it must 
be noted that Lagrange coordinate q in the above argument is not the initial position of 
the sheet. Please refer to Roytvarf (1987), and Gouda & Nakamura (1989), for the detailed 
explanation. On the other hand, it is found that p(x) is proportional to x~ 2 ^ 3 around the 
singular points at the first appearance of caustics (Zel'dovich (1970), Arnold et.al.(1982), 
Gouda & Nakamura (1988)). This type of singularity is called A3 type. The A3 Type is 
not structural stable in the one-dimensional system and then the A3 type of singularities 
appears only at an instant. It disappears and quickly evolve into the A2 type, i.e. p oc x^ 1 ^ 2 
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singularity. Hence the contribution from the A3 type singularity in estimating the power 
spectrum is negligible. The fourier spectrum of the density is given by 

S k = J 5{x)e ikx dx 

^e lkx dx 
Po 

Ikx , 



where t = kx. Then, we obtain 



oc J (f3x)—ze lkx dx 
= p-hh' 1 Jt-*e*dt, (12) 



P{k) oc (k(3)-\ (13) 

From eq.flilf), we find the value of the power index is —1. Here we notice that small (3 
contributes to the large amplitude of the power spectrum. Furthermore, the coefficient (3 
is very small when the caustics have just appeared. We can see this fact by the evolution 
of the "single-wave perturbation" . Here, the single-wave perturbation means the density 
fluctuation whose initial condition is given by 

x (l) — 1 + Xsin(q), (14) 

where X is a constant value. Hereafter we call this case the single-wave case. The first 
caustics have just appeared at a center of the x axis (x = 0) in Fig. 4(a). After the first 
appearance of this caustic, more caustics appear in the collapse regime by phase mixing 
(Fig.4(b)-(d)) (Doroshkevich et.al(1980), Melott (1983), Gouda & Nakamura (1989)). A 
caustic (singularity of the density field) is located at the point where the derivative of v 
with respect to x is infinity. The absolute value of the derivative of v with respect x around 
the singular points is proportional to (3~ l . Then as the phase mixing continues and so 
more caustics appear, it is found from Fig.4(a)-4(d) that (3 increase because |^| around the 
singular points decrease. Hereafter we call the structure in the phase space shown in Fig.4 
(a)-(b) "the whirlpool". Just after the multi-caustics have appeared, the coefficient factor 
P of this wavelength is very small and the effect of the amplitude of the power spectrum in 
these waves dominates. As a result, the index of —1 that is predicted by the catastrophe 
theory appears. Therefore we call these regimes, single-caustic regime. 

Here, we comment that in the real singularities do not occur in the real world. Of 
course, the arguments on the caustics in this paper cannot be applied to barionic gas 
because shocks appeared and they prohibited that the density diverges. Furthermore the 
maximum density is surprisingly low even in the neutrino-dominated model due to their 
velocity dispersion (Zel'dovich & Shandarin (1982),Kotok & Shandarin (1987)). In this 
paper, we consider the cold collisionless matters whose velocity dispersion is much small 
compared with the scales under consideration. 
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2.2.2. Multi- caustics regime (the regime 3) 

In the multi-caustics regime (k > k sn i), we can see the power law spectrum. With 
the power index of these regimes is different from —1 predicted by catastrophe theory. In 
this regime, various small scale fluctuations (k ^> k sn i) have already collapsed and made 
singularities. Every singularity makes "the whirlpool" in phase space as shown in the 
single- wave case (Fig. 4). Various size of "whirlpools" are made of the various scale of 
fluctuations. The distribution of the whirlpools determine the value of the power index. 
This distribution depends on the initial power index n. That is, the power index /i in 
this regime depends on n. Indeed, the value \i depends on the initial power spectrum as 
shown in Fig. 2 (a), and 2(b). Therefore we call this regime the multi-caustics regime. When 
the self-similarity is satisfied, there is a relation between the power index of the two-point 
correlation function, the initial power spectrum index, and the mean relative peculiar 
velocity. Here, we briefly show this relation in the w dimensional system. In the linear 
regime of the w dimensional system, the two-point correlation function is given by(Peebles 
1980,1993), 

f oc a 2 x-^ oc (^)- {w+n \ (15) 

where 

a = . (16) 

w + n 

On the other hand, in the strongly nonlinear regime, the two-point correlation function 
obeys the following evolution equation (0th moment equation of the 2nd BBGKY eq.: Davis 
& Peebles (1977), Yano & Gouda (1997)), 

§ + ^T^ 1, « + 1 ><^ ' < 17 > 

In the strongly non-linear limit, we can assume a power law solution for the two-point 
correlation function (Yano & Gouda (1997)). 



£ oc a e x~^. (18) 

In this case, we obtain 

( v ) = -hax, (e - h{w - 7) = 0) (19) 

and then, 

$ oc a hiw ^x^ = (JL)-t, (20) 
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where h is the constant and 

a' es KW - 1] . (21) 

7 

Hereafter we call the parameter h the relative velocity parameter. If the self-similarity is 
satisfied, a' must be equal to a. Then the power index 7 is given by 

' 2 + h{n + w) K J 

If the relative velocity parameter h is 1 (stability condition) as Davis & Peebles (1977) 
assumed, we obtain the 7 = (n + l)/(3 + n) in the one-dimension (w = 1). In this case, the 
power index \x of the power spectrum is given by // = 7 — w = —2/(3 + n). However the 
power index fi which is given by the numerical results of the power spectrum is different 
from — 2/(3 + n). Therefore, the stability condition (h = 1) is not satisfied. We can estimate 
the velocity parameter h from the power index by using eq.( 



h 2 7 + 2Q/+1) 

(n + w)(w — 7) (n + w)fi (n + l)fi 

We show the results in the Table 1. Yano & Gouda (1997) discussed that the velocity 
parameter h has the value between and 1. Indeed the value of h stays at this range. We 
show the numerical result of h not only for the n = 1, and 2 cases but also for the n = 0, 
and 3 cases. 

Table 1. n dependence of \x and h 

n fj, h 

^88 029 

1 -0.75 0.33 

2 -0.65 0.36 

3 -0.58 0.37 



2.3. Cutoff-case 

We consider the following initial spectrum, that is, the cutoff-case. 

^»{^o Ifdr (24) 

The time evolution of the power spectrum for the cutoff-case is shown in Fig. 5. It is shown 
only for the n = 1 case. However, we obtain essentially the same result for other cases 



-10- 



(n = 2,3, and 0). As we can see from Fig.5, until the first appearance of caustics (lowest 
solid line of the power spectrum except for initial spectrum in Fig.5), the amplitude of the 
power spectrum grows even on scales smaller than the cutoff scale. The power index of 
these regime is predicted to be — 1. Because the wave number k = k cu t makes the caustic 
and this caustic results in the predicted value of power index as shown in §2.2. Indeed, we 
find that the value of the index in these regimes is nearly equal to —1. The scaled power 
spectrum is also shown in Fig. 6. The definition of the scaling is the same as that in the 
scale-free case. As we see from Fig.6, each power spectrum at the different time does not 
coincide. Therefore the self-similarity is not satisfied. 

The characteristic scale of separation of the caustics becomes smaller and smaller as 
time increases (refer to the case of the single- wave: See Fig.4 (a)-(d)). Here, k cs is defined 
as the wave number of the characteristic scale. On the scales smaller than the characteristic 
scale k > k cs (regime 5: See Fig. 10), the power index of the power spectrum becomes —1. 
Because on the scales smaller than the characteristic scale, the density profile around one 
singular point determines the power spectrum. Therefore we call this regime the smallest 
single-caustic regime. On the other hand, the regime k cu t < k < k cs (regime 4), the 
power index of the power spectrum is different from — 1. The index has a certain value v. 
Because on the scales larger than the characteristic scale, smoothed density profile with the 
smoothing scale {k cut < k < k cs ) determine the power spectrum. We call this regime the 
virialized regime. We show the schematic of the density distribution in Fig. 7. We notice 
that the wave number k cs which represents the characteristic separation of caustics becomes 
larger and larger as time increases. On the contrary, the non-linear scale wave number k n i 
and k sn i become smaller and smaller as time increases. Thus, the power spectra of the 
different times cannot coincide on all scales. Therefore the self-similar evolution in the all 
regimes including the scales smaller than the cutoff scale cannot be satisfied. 

We compare the time evolutions of the power spectrum for the cutoff-case with those 
for the scale-free case in Fig. 8. Of course, the power spectrum for these two cases does not 
coincide with each other on all scales. On the other hand, on the scales larger than the 
cutoff scale, we can see the coincidence of power spectra for these two cases. This means 
that the growth rate on the wave number smaller than k cu t does not affect the difference of 
the larger wave number than the cutoff scale. We show the scaled power spectrum for the 
cutoff-case only on the wave numbers smaller than the cutoff scale k cut in Fig.9. This means 
that even in the cutoff-case, the self-similar evolution can be satisfied on scales larger than 
the cutoff scale. 
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3. CONCLUSIONS AND DISCUSSION 

We have calculated the time evolution of the power spectrum for two cases of initial 
conditions. One is the scale-free case and the other one is the cutoff-case. In the case of 
the scale-free case, we can see the self-similar evolution of the power spectrum. The scaled 
power spectrum P*(fc*) at each time coincides with each other. We can separate roughly 
the power spectrum into three regimes. One is the linear regime (k < k n \ : the regime 1). 
The value of the power index in this regime remains the initial power index, n. The second 
regime is the single-caustic regime (k n i < k < k sn i : the regime 2). In this regime, the 
power index becomes —1 and is independent of the initial conditions. This result is caused 
by the appearance of caustics at this scale, and these caustics determine the power index of 
the power spectrum in this regime. The third regime is the multi-caustics regime (k > k sn i 
: the regime 3). The distribution of the "whirlpool" in phase space determine the value of 
the power index. Therefore in this scale the power index \i has the value which depends on 
the initial condition. We can estimate the velocity parameter h from the power index in the 
multi-caustics regime, which is around 0.5. Yano & Gouda (1997) discussed the probable 
value of h from the physical point of view, and obtained that h takes a value between and 
1. Estimated values of h are consistent with this argument. Indeed, the index /x and the 
velocity parameter h depend on the initial power index n. The stability condition (h = 1) 
is not satisfied in this case. 

In the cutoff-case, we find that there is no self-similarity on all scales. The scaled 
power spectrum does not coincide with each other in all regimes. However, the spectrum 
coincides on the scales larger than the cutoff scale. After the appearance of the first 
caustics, the power index in the regime of the scale smaller than the cutoff scale becomes 
— 1. This value is, as mentioned above, caused by the appearance of caustics as shown in 
Fig. 3. More and more caustics appear one after another and so the separation of caustics 
becomes smaller and smaller. Even after the appearance of many caustics, on scales smaller 
than the characteristic separation of the caustics (k > k cs : the regime 5), the power 
index of the power spectrum is obtained by the density profile around the singular point. 
Therefore the power index of these scales becomes —1 which can be derived according to 
the catastrophe theory. On the other hand, on the scales larger than the characteristic 
separation of the caustics (k cut < k < k cs : the regime 4), the distribution of the singular 
points determines the power index on these scales instead of the density profile around 
one singularity. Because in this regime, the smoothed density profile with the smoothing 
scale (k cut < k < k cs ) determine the power index. On these scales, the distribution of the 
singularity occurring in the evolution of the single- wave is important. Then, the power 
index on these scales is determined by this distribution of the singularity. Therefore, we 
can roughly separate two regimes on the scales smaller than the cutoff scale k cut . One is 
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the virialized regime (k cut < k < k cs : the regime 4), and the power index has the value v . 
Another is the smallest single- caustic regime (k > k cs : the regime 5) and the index have 
the value of —1. Both indexes are independent of the initial conditions. 

In a real situation about the evolution of the power spectrum, the initial power 
spectrum has a cutoff at a certain scale. Therefore, in the case of a real situation, we 
can consider the evolution of the cutoff-case. We notice that even if the initial power 
spectrum does not have a cutoff, it decreases with the power index n < — 1 on certain scales 
(k > kdec)- Therefore the time evolution of the power spectrum is the same with the one in 
the cutoff-case after the appearance of caustics. And then, we can separate roughly five 
regimes. We show a schematic of the power spectrum at a time after the first appearance 
of caustics in Fig. 10. First regime is the linear regime (k < k n i : the regime 1). The index 
of the power spectrum in this regime is n. Second regime is the single- caustic regime 
(kni < k sn i : the regime 2). The index in this regime is —1, and independent of the initial 
conditions. Third regime is the multi-caustics regime (k sn i < k cut : the regime 3). The 
power index is \i which depends on the initial power index n. Fourth regime is virialized 
regime (k cut < k cs : the regime 4). The value of the power index is v which is independent 
of the initial power index. The fifth regime is the smallest single-caustic regime (k > k cs : 
the regime 5). The value of the power index is —1. Only when we consider the evolution of 
the power spectrum on scales larger than the cutoff scale (i.e. the first, second, and third 
regime), the self-similarity is satisfied. 
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Fig.l.-(a) The time evolution of the power spectrum. The solid and dotted lines are used 
mutually in order to distinguish easily each line. The curves shown are (bottom to top) 
at a = 0.1, 1, 2, 4, 8, 16, 32, 64, 128. The initial power spectrum is scale-free with the power 
index n = l(Dotted-dashed line). 

(b) The same as (a) but for the initial power index n = 2 and the power spectra are shown 
at the scale factor a = 0.01, 1, 2, 4, 8, 16, 32, 64, 128. 

Fig. 2. -(a) The scaled power spectrum P*(fc*). The initial power spectrum is scale-free with 
the power index n — 1. Three solid straight line shows the power law with the power index 
of 1,-1,-0.75. 

(b) The same as (a) but for the initial power index n — 2. Three solid straight line shows 
the power law with the power index of 2,-1,-0.65. 

Fig. 3.- The schematic of the power spectrum in the scale-free case at a certain time. 

Fig. 4.- (a) The distribution in the phase space for the single- wave case at the first appearance 
of caustics. A horizontal axis and a vertical axis are the Eulerian coordinate x and the 
velocity, respectively. Scale factor a is normalized at this time (a = 1). 

(b) The same as (a) but at a = 4. 

(c) The same as (a) but at a = 16. 

(d) The same as (a) but at a = 32. 

Fig. 5.- The time evolution of the power spectrum. The solid and dotted lines are used 
mutually in order to distinguish easily each line. The curves shown are (bottom to top) at 
a = 0.05, 1, 2, 4, 8, 16, 32, 64, 128. The initial power spectrum obeys the power-law with the 
index n = 1 and it has a cutoff at k cut = 127(Dotted-dashed line). 

Fig. 6.- The scaled power spectrum P*(k*). The initial power spectrum obeys the power-law 
with the index n = 1 and it has a cutoff scale at k cut = 127. 

Fig. 7.- The schematic of the density distribution in the single- wave case. This figure shows 
that the density profile around one caustic and the distribution of caustics, which determine 
the smoothed density profile on scales larger than the characteristic separation of caustics. 

Fig. 8.- Comparison of the evolution of the power spectrum between the scale-free case (solid 
line) and the cutoff-case (dashed line). Both the solid curves and dashed curves shown are 
(bottom to top) at a = 0.05, 1, 2, 4, 8, 16, 32, 64. The initial power spectrum obeys the power 
law with the power index n — 1 both for the scale-free case and for the cutoff-case. For the 
cutoff-case, the initial power spectrum has a cutoff at k cut = 127. 
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Fig. 9.- The scaled power spectrum at a=l, 2,4, 8, 16, 32, 64 on scales larger than the cutoff 
scale. The initial power spectrum obeys the power law with the power index n=l and it 
has a cutoff scale at k cut = 127. 

Fig. 10.- The schematic of the power spectrum with a general initial power spectrum at a 
certain time after the first appearance of caustics. 
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